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BARYON LOADING OF AGN JETS MEDIATED BY NEUTRONS 
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ABSTRACT 

Plasmas of geometrically thick, black hole (BH) accretion flows in active galactic nuclei (AGNs) are gen- 
erally collisionless for protons, and involve magnetic field turbulence. Under such conditions a fraction of 
protons can be accelerated stochastically and create relativistic neutrons via nuclear collisions. These neutrons 
can freely escape from the accretion flow and decay into protons in dilute polar region above the rotating BH 
to form relativistic jets. We calculate geometric efficiencies of the neutron energy and mass injections into the 
polar region, and show that this process can deposit luminosity as high as Lj ~ 2 x 1(T 3 Mc 2 and mass loading 
Mj ~ 6 x 10~ 4 M for the case of the BH mass M ~ 10 s M Q , where M is mass accretion rate. The terminal 
Lorentz factors of the jets are Y ~ 3, and they may explain the AGN jets having low luminosities. For higher 
luminosity jets, which can be produced by additional energy inputs such as Poynting flux, the neutron decay 
still can be a dominant mass loading process, leading to e.g., r ~ 50 for L iitot ~ 3 x 10~ 2 Mc 2 . 
Subject headings: galaxies: jets — black hole physics — plasmas 



1. INTRODUCTION 

One of the major problems in astrophysics is the produc- 
tion mechanism of relativistic jets. They are associated with 
active galactic nuclei (AGNs), Galactic black hole (BH) can- 
didates, and gamma-ray bursts (GRBs). It is inferred that 
AGN jets have Lorentz factors of Y ~ 10- 100, luminosities 
as high as the Eddington luminosity Lndd, and opening angles 
of 6j ~ T" 1 . Although the matter content of AGN jets is still 
an op en problem, the in ertia is seemingly dominated by pro- 
tons (ISikoraet al.ll2005l) . The masses of central BHs of AGNs 
are typically M ~ 10 7 M Q - 10 9 M Q . 

An outflow can have a relativistic velocity if the enthalpy 
per unit rest energy /i = Lj^jMp 1 ^ 1 in the vicinity of the 
BH, where L /tot = Lg+Lk is the total luminosity, Lb and 
are Poynting and particle kinetic luminosities, respectively, 
and Mj is mass loading rate. If the source of L 

j, tot is gravi- 
tational energy of the accreting mass on the central BH, we 
have L / tot < GMM/R < Mc 2 , where M is mass accretion 
rate, and radius R should be larger than GM/c 2 . This leads 
to M /Mj 3> 1, implying that the production mechanism of 
relativistic jets should be converting gravitational energy into 
Poynting and/or kinetic energies, and concentrating them on 
a small fraction of mass. 

It is therefore likely that the Poynting and/or kinetic en- 
ergies are injected in polar region above the rotating BH, 
'the funnel', where mass lo ading is exponentially sup pressed 
by the centrifugal barrier (Abra mowicz et all 11978b . This 
configuration is also suitable for the outflow to be colli- 
mated by the external pressure. The magnetically dominated 
jet models, or magnetohydrodynamic (MHD) models, have 
been recently pr ogressed with numerical simulations (e.g., 
lMcKinne\l2006l) . If a certain amount of particles are injected 
(or generated) in the funnel, the strong poloidal magnetic 
fields associated with electric currents flowing in the accretion 
torus accelerate the flow of the particles. The MHD flow can 
accelerate to a relati vistic velocity if suitable boundary condi- 
tions are satisfied dKomissarov et all 120071; lLvubarskvll2009t 
iGranot et al.ll2.01 lb . On the other hand, the kinetically domi- 



nated jet models, so-called fireball models, have been actively 
discussed for GRB jets, for which the energy injection via 
vv an nihilation could be efficient (Zalamea & Beloborodov 
2011). A thermally dominated spherical blob of gas in- 
evitably accelerates to a relativistic velocity (Meszar os et al.l 
fl993tlPiran et al.H993t lkobavash i et alJ1999ft . It is even pos- 
sible that large amount of remaining thermal energy at the 
photosphere is re leased as prompt 7-ray emission itself (e.g., 
iToma et aill201 ll and references therein). For AGN jets also, 
firebal l models have been elaborately studied by some a uthors 
(e.g.. lAsaiio & Takaharall2007l 120091; [Becker et al.ll20Tll) . 

However, the above models do not answer a question "Why 
are the Lorentz factors of AGN jets regulated to Y ~ 10- 100? 
Why not Y » 100 or Y ~ a few?" In either of the MHD 
or fireball model, mass loading rate into the funnel is essen- 
tial for determining the final Lorentz factor (and the radia- 
tion properties) of jets. There is an interesting idea on this 
point for GRB jets. Their mass loading might be determined 
by neutron diffusion from the disk wind sur rounding the jet 
dLevinson & Eichlerll200l lMcKinnevll2()05T) . The jet is con- 
sidered to have globally ordered magnetic fields that may sup- 
press proton diffusion across the jet-disk wind boundary^ 

In this paper, we discuss the role of the neutrons for 
mass (as well as energy) injection of AGN jets. We fo- 
cus on relativistic neutrons escaping from the accretion flow, 
and calculate the fraction of the total escaping neutrons 
that decay into protons in the funnel. The processes of 
the relativistic neutron escape from the B H accretion flow 
have been discussed by several authors (e.g..lEichler & Wiital 
[19781; ISikora et dJH989t rBegelmanetalJl 9901; lAtovanll 19921; 
iContopoulos & Kazanasn l995). although the efficiencies of 
mass and energy injections in the funnel region have not been 
studied in detail. In Section [2] we make a short review on 
proton acceleration and cooling in the accretion flow, and 
parametrize the spectrum of the escaping neutrons. Then we 
formulate the geometric injection efficiency in Section|3] and 
show the results of Lj/Mc 2 and Mj/M in our model in Sec- 



1 See also Levinson ( 2006) and Ioka 12010) for other ideas on relativistic 
jet mass loading. 
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tionlU Conclusion and discussion are given in Section[5] 

2. RELATIVISTIC NEUTRON PRODUCTION AND ESCAPE 

We consider the AGN central engine as follows. A rotat- 
ing BH has an accretion flow around it. The accretion flow 
is geometrically thick, and creates a hot corona and/or non- 
relativistic wind extending vertically, although the polar re- 
gion above the central BH is dilute because of the centrifugal 
barrier. There can be large-scale magnetic fields anchored on 
the accretion flow, which may also suppress protons leaking 
into the funnel, since the Larmor radius 



eB 



3.1 x lO^pBJ 1 cm 



(1) 



is much smaller than the system size characterized by the 
Schwarzschild radius 



R s = 



2GM 



3.0 x 10 1J M 8 cm, 



(2) 



where B 3 = B/(1Q 3 G) and M 8 = M/(10 8 M Q ). Neutrons are 
not affected by the magnetic fields, so that they can be an 
effective source of mass injection into the funnel. 

One may consider the origin of such neutrons as thermal 
process in the accretion flow. The temperature of protons and 
heavy nuclei is as high as kT p ~ GMm p /R ~ 500/?; 1 MeV, 
where = R/R s . Then the helium breakup and the pp col- 
lisions can produce thermal neutron component. The neu- 
tron fracti on can be n„/(n n + n n ) 5, 01 in th e vicinity of the 
BH (e.g.. lFilho et~aT]l2003h iHu & Pengll2008l) . However, non- 
relativistic neutrons only survive as long as 



D = ct„ ~3 x 10 13 cm, 



(3) 



where r„ ~ 10 3 s is the decay time of a neutron. This is com- 
parable to the system size characterized by R s , and thus they 
cannot reach the funnel in the case of M > 10 8 M Q . 

We then focus on relativistic neutrons created in the ac- 
cretion flow. The plasmas of the geometrically thick AGN 
accretion flows are generall y collisionless for protons (cf. 
Takahara & Kusunose 1985). In such plasmas, the pro- 
ton energy distribution is not necessarily Maxwellian, and 
some fraction of the protons can be accelerated to relativis- 
tic speeds. The relativistic hadrons produce relat ivistic neu- 
tron component via the pp a nd/or 777 collisions ( Siko ra et ail 
Il989t Begelma n et al.lfl990h . The relativistic neutrons with 
Lorentz factor 7„ survive over D^f„ ~ 3 x 10 13 7„ cm, so that 
they can reach the funnel and decay into relativistic protons 
(and electrons) there. The magnetic coupling in the funnel 
will isotropize them, leading to the electron-proton fireballs. 
Neutrons that decay outside the funnel just energize the dense 
non-relativistic wind. 

As for the proton acceleration mechanism in the geometri- 
cally thick accretion flows, there are many possibilities, and 
we do not specify it in this paper. If the accretion results from 
turbulent viscosity drive n by magneto-rotational instability 
(Balbus & Hawlev 1991), the plasma involves magnetic field 
fluctuations and reconnections, where protons can be acceler- 
ated via second-orde r and/or first-ord e r Fermi processes (e .g ., 
iDermer et aD 1 19961 iHoshinol 120121: iRiquelme etai] 120121) . 
The flows may also have regions with strong (less-fluctuated) 
magnetic fields, where accretion can result from angular mo- 
mentum transport by magnetically driven winds. Some fields 
may be oppositely directed, giving rise to magnetic recon- 
nections, which can accelerate particles via first-order Fermi 



process (e.g.. lDrurvll2012l Ide Gouveia Dal Pino et all 120101: 
IVievro & Romeroll2012l) . Furthermore, protons could be ac- 
celerated at a stand ing shock formed in the accretion flow 
(lBeckeretal.1120111) . 

In th e f ollowing, we make a sh ort review of Siko ra et al.l 
( 1989) and Begelma n et al.l {1990) for parametrizing the neu- 
tronization factor and the escaping neutron spectrum. 

2.1. Spectrum of escaping neutrons 

Relevant processes for protons are the accretion inflow, 
Fermi acceleration, pj cooling, and pp cooling, and those for 
created neutrons are the escape from the inflow, n-f cooling, 
and np cooling. The proton escape or spatial diffusion is neg- 
ligible because of the small Larmor radius, whereas protons 
are advected on the inflow timescale fj n : 



R , R 



(4) 



where v; n = 0.1v; n ,-ic is the accretion fluid velocity. Fermi 
acceleration timescale is estimated as f acc ~ £Rl/c — 1.0 x 
10~ 3 jp^B^ 1 s, where £ = 10 4 ^4 is the factor related to the 
magnetic field fluctuations. This indicates that the accelera- 
tion is so efficient that protons can be highly energized. At 
high energy ranges, the pj collisions are the most important 
cooling process, where we assume a typical radiation field 
with luminosity 10~ 4 LEdd < L r < Lgdd and spectrum F v oc v~ a 
with 0.5 < a < 1.5. This assumed radiation field is rather 
generic, which may include thermal as well as non-thermal 
emission of electrons. The equality of f acc with the p^ cool- 
ing timescale t pl leads to the maximum proton Lorentz factor 

lpM ~ 10 7 (B 3 /?*/e4) 1/2 (i,/10- 2 L Ed dr 1/2 . 

For the BH mass M ~ 10 7 M© - 10 9 M©, the system size 
is < 10 15 cm, so that the protons with j p < 10 2 (which are 
converted into neutrons with similar Lorenz factor 7,,) are im- 
portant. At such energy ranges, the pp collisions are more 
efficient for creation of neutrons than the 777 collisions. The 
pp cooling time for protons is estimated as 



1 



PP v- 

n p c(7 pp J\ pp 



' 30 rZ 



,R 



(5) 



where n p and t p = ojn p R are the proton number density and 
the Thomson optical depth at radius R, respectively. The pp 
cross section is a pp ~ oj j 17, andK pp ~ 1/2 is the inelasticity. 
If t pp < f; n , the pp collisions occur efficiently. For t m < t pp , 
the efficiency reduces by a factor of t- m /t pp . 

The created neutrons are not magnetically coupled to the 
background plasma. They escape without being absorbed if 
f e"c ~ ^/c is smaller than the timescales of np and nj col- 
lisions, t„ p and f„ 7 , that are similar to t pp and t pl , respec- 
tively. The maximum Lorentz factor of the escaping neutrons 
is given by 4c = h n , as j nM - lO 6 ^/^" 2 ^)" 1 . 

At lower energy ranges, say j p < 10 5 , we have f acc < < 
min(ti n ,t pp ) < max(ti n ,tp P ) < t PJ , for t p < 20. In this case, 
neutrons created by pp collisions freely escape from the ac- 
cretion flow. The kinetic equations for the proton and neutron 
number densities in the steady state and for a power-law pro- 
ton injection function N p oc j~ p indicate that the ratio of the 

escaping neutron number flux ^V„ to the proton injection flux 
is estimated as 



1 



N p 2^-D + l 



1 1 ln 
mm 1 , — 



'/'/' 



(6) 
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The numerical factor 2 in the denominator is determined by 
the probability of the charge exchange during a single pp col- 
lision and K pp . The neutron spectrum is given by N n oc -f~ p 
for the low-energy range satisfying t pp < t pl . It deviates from 
a single power-law at hi gh energy, say j„ > 10 5 (see Fig. 7 
of Begelm an et al.l [l990). although the neutrons at such high 
energy ranges are not relevant for our purpose in this paper. 

The process of the creation and escape of neutrons with 
7„ < 10 5 is most effective for < t np and t pp < f in , i.e., 

3v in _! < t p < 20, (7) 

for which we have /„ ~ 1 /[2{p- 1)+ 1]. It may be possible 
that geometrically thick accretion flows satisfy this condition 
of T p , To confirm it in detail, however, modeling of accre- 
tion flows with significant neutron energy release is required, 
which we leave as separate work. 

3. CALCULATION OF THE GEOMETRIC EFFICIENCY 

The luminosity of the escaping neutrons is parametrized by 
L„ = f„f a faMc 2 , where fa is the ratio of the heating rate of 
protons to Mc 2 , f a is the ratio of the rate for accelerated pro- 
tons to the heating rate of protons, and /„ is the neutronization 
factor estimated by Eq. ©. The heating rate can be as high 
as - (l/2)GMM/R s ~ Mc 2 /4, i.e., fa < 0.3, while we have 
fn ;$ 1 /[2(p— 1)+ 1] < 0.3 for a reasonable range p > 2. It 
is difficult to estimate the acceleration fraction f a . We only 
have a rough constraint f a < 0.3, which means that the en- 
ergy density of the accelerated protons will not dominate that 
of thermal protons. Therefore we may summarize the micro- 
physical efficiency as 

Ufa fa < 3 X lO- 2 . (8) 

The neutrons with this luminosity are released isotropically, 
and the protons created via the neutron decays are magneti- 
cally coupled and energize the background plasma. Here we 
calculate the geometric efficiency, i.e., the fraction of neutrons 
that decay into protons at the polar region. In order to obtain 
the order of magnitude of the geometric efficiency, we set a 
simple configuration of the BH accretion system, and calcu- 
lations are performed by assuming the Euclidean space. 

As illustrated in Figure Q] we consider an accretion torus 
with H/R=l/ tan 6>o around a BH, which is located at the cen- 
ter of the coordinate system. The neutron emission region is 
assumed as the inner fraction of this torus between the cylin- 
drical radii p m = and pm = 3R S . We assume that the neutron 
emissivity is uniform in this region for simplicity. The spatial 
volume of this region is V = (4/3)7r(p^- p^/tan^o- Then 
the neutron emissivity is given by dE ( "^ /dtdVdQdj = ffi = 

A 7 "/> +1 for 1 < 7 < 7 M, where A = L„/(4ttV J^' -f~P +1 d-f). 
Hereafter we will not distinguish the Lorentz factors of pro- 
tons and neutrons, since the neutrons with j„ decay into pro- 
tons with Tp ~ 7„. The dilute polar region is set as the cone 
with the opening angle 9\ . We assume that 0\ ~ 0(1) for 
n = R s < Tj < r%, and consider that r2 is a characteristic ra- 
dius above which the collimation becomes strong due to the 
pressure from the hot corona and/or disk wind. The collima- 
tion suppresses the energy injection for rj > ro- We will take 
r2 = 10R S as a fiducial case below. 

Let us consider a straight line with a parameter s cross- 
ing the surfaces of the emission region at s = s\ and s = s>i. 
The fractional neutron intensity emitted from the line ele- 
ment As at si < s < S2 and measured at a point s in the 




FIG. 1. — Structure of the BH accretion system for the calculation of the 
geometric efficiency. The BH is at the center of the coordinate system. The 
accretion flow with thickness H/R = l/ tan 8q has the neutron emission region 
between the cylindrical radii p„, = R s and pu = 3R S . A trajectory of neutrons 
incident to the point r, = (rj,6j,ipj) is shown. The dilute polar region is 
assumed to be a cone with the opening angle 9\ for rj < ("2 = 1(M? S , and 
strongly collimated for rj > rj. 

polar region is Al^(s) = f"\s)As, if neutron decay is not 
taken into account. We may approximate the proton inten- 
sity created at Sj < s < sj + dsj through the neutron decay 
as dAl^/dsj = Al^\s)5(sj-s-Dj). Integrating this over 
Sj < s< S2, we obtain the proton emissivity created at Sj from 
the neutrons propagating through a single line as 

qj(p) r*2 
j^(sj) = -^(sj) = J f;\s)6(Sj-s-D 7 )ds 

= f A-f~P +l (sj -s 2 <Df<Sj-s{), (g) 
[0 (otherwise) 

For a given point r = rj (corresponding to Sj), S\ and S2 are 
functions of the incident direction, i.e., s t = Si(f2j) and «2 = 
S2(£lj). We obtain the energy injection rate of protons per 
unit volume at a point r = r ; by integrating ffl(Sj) over solid 
angle, dE^ /dtdVjdj = u^\rj) = J dQjj^\sj). This leads to 
the mass injection rate per unit volume as dM^ p) j dtdVjd'y = 
w'^(r ; )/(7c 2 ). Finally we obtain the total energy and mass 
injection rates by integrating and «^ ) /( 7 c 2 ), respectively, 
over the total energy range and the volume of the polar region . 
The efficiencies of the neutron energy and mass injections are 
thus written by 

rv(M C 2 )i_^ f ei de . sin0 . r dr r 2 X 

\Mj/M )-Mc 2 L '"'A. ' ;X 

I d V da '\khY am 

where the axisymmetry of this system has allowed us to per- 
form the integration over the azimuthal angle ipj, and we set 
the upper bound of the rj integration as r2 since the energy 
injection is assumed to be much less effective for r, > r2 due 
to strong collimation. The free parameters for calculating the 
efficiencies for a given BH mass M are tan0OiPi7M, r 2,tan#i, 
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lj=1 lj=3 rj=10 rj=30 rj=60 




1 10 100 

7 

FIG. 2. — Fraction of spectral energy injection rate per unit volume Wy(rj) 
(times rj) of protons at the pole (0j = 0) of various radii rj, with respect 

to the escaping neutron luminosity L„ = fnfafthMc 2 , for the case of Mg = 1, 
The lines from left to right correspond to rj = 1,3, 10,30, and 60, respectively 
(measured in unit of D). The other parameters are tan 6q = 2,p = 2, and -yu = 
10 2 . 




1 10 100 

Y 



FIG. 3. — Fraction of spectral energy injection rate per unit solid angle for 
1< rj < r 2 = 10 dLj/d(cos6j)d<p j d'y = J[ 2 drj rjiif } at the pole (0, = 0), 
with respect to L„ = fnfafthMc 2 , for the case of Mg = 1 with tan#o = 2 and 
7M = 10 2 . The solid, dashed, and dot-dashed lines correspond to the cases 
of p = 2, p = 3, and p = 1, respectively. For reference, we plot a power-law 
function oc 7 _1 by the solid straight line. 

and fnfafth- 

4. RESULTS 

4.1. Case of M= 10 s M 

Here we show the calculation results for the case of Ms = 1 . 
In this section we measure lengths in unit of D. Since R s 
happens to be equal to D in this case, we have p m = 1 and 
Pm = 3. The other parameters are set to be tan#o = 2,p = 
2,7m = 10 2 ,r2 = 10, and tan 9\ = 1. We will examine the de- 
pendence of results on these parameters later. First, we show 
the integration result of ifflirf) = j diljj^Xsj) to see the spec- 
tral property of the injection rate at various points. We plot 
rjii^' /L„ at the pole 9j = of various radii rj = 1,3, 10,30, 
and 60 in Figure [2] We can see that the injected protons are 
distributed over relatively broad energy range at r, < 3, while 
concentrated to 7 ~ r, at large radii rj > 3. This indicates that 
for large r,, neutrons decay after traversing large distances, so 
that the detailed structure of the emission region is not rele- 
vant. 



0.05 r 

0.04 - 



«- 0.03 - 




jjf 0.02 - 
_j 

0.01 - 

I ' • ' • ' • ' 1 ' 

0.2 0.4 0.6 0.8 1 

8j [rad] 

FIG. 4. — Fraction of energy injection rate per unit solid angle for 
1< rj < r 2 = 10 (dLj/d(cos6j)d(pj = dy f[ 2 drj rju^) with respect 
to L„ = fnfafthMc 2 as a function of 8j for the case of Mg = 1 with the other 
parameters tan 9q = 2,p = 2, and 7m = 10 2 . 

Next we show the integration result of J^ 2 drj rju^ at the 
pole, which is equivalent to dLj / d(cos9 j)d(p jdj . We plot this 
value, normalized by L„, in Figure [3] We also show the re- 
sults for p = 1 and p = 3. For p = 2, comparing the result 
with the injection power-law profile oc j~ p+l , we find that the 
energy is injected efficiently for 7 ~ 3 — 10, and has the simi- 
lar spectrum as the injected one. This property is explained 
simply as follows. At the large radii, say rj > f\, where 
the detailed structure of the emission region is not relevant, 
the solid angles of the incident neutron directions are limited 
within a small range A£lj ~ S/rj, where S is the horizon- 
tal cross section of the emission region. Also we may write 
j^'irj) = j^'8{rj — r i)h, where h represents the mean width of 
the emission region. Note that we measure r, in unit of D. 
These approximations lead to 

— f dr, r 1 :U { !' ) ~ -tJ- x 

L n J h 1 1 7 wtf M r p+l di 

J drjPj J dnjj- p+] S(rj-j)h 

r^v^~' ,+I (?i<7<r2) ai) 

1.0 (otherwise), 

where we set Sh ~ V. Since we have 1/(4tt J™ 7-^7) = 

0.017 for p = 2 and jm = 10 2 , this rough calculation well 
agrees with the numerical integration for 3 < 7 < 10. For 
7 < 3, the spatial volume of emitting points that can connect 
to the pole with decay lengths of ~ 7 is smaller, so that the 
injected energy is smaller. 
We integrate j- f{ 2 drj rju^ over 1 < 7 < j M = 10 2 and 

then obtain j-dLj/d(cos9j)dipj = 0.023 at 9j = 0. We also 
show the calculation results of j-dLj/d(cos9j)dtpj as a func- 
tion of 9j in Figure|4] This shows that the energy injection rate 
does not strongly depend on 9j. The slight increase for larger 
9j is due to larger contribution for the energy range 7 < 3, but 
this contribution is confirmed to be not significant for the total 
energy injection Lj in the case of tan(9i = 1. 

Finally we integrate j-dLj /d(cos9 j)dipj overO < 9j <9\ = 
7r/4 and multiplied by 2tt, obtaining Lj/L„ = 0.051. This 
provides the total energy injection efficiency for the param- 



Relativistic neutrons for AGN jet production 



0.25 - 



0.15 




0.05 - 



FIG. 5. — Factor J* 3 '°7 p+1 ^7/ J^ M 7 ,,+1 rf7 is plotted as a function of p 
for various 7« . 

eter values tan(9o = 2,p = 2,7m = 10 2 ,r2 = 10,tanf?i = 1 as 
Lj/(Mc 2 ) = 0.051 f„f a f th . 

The geometric efficiency is roughly estimated by integrat- 
ing Eq. (fTTT> with c/(cos 9j)dipjdj, 



Lj_ 2tt(1-cos6>i) J^7 p+l d-f 



(12) 



which gives ~ 0.04, if ri = 3 is taken. The slight difference 
between this rough estimate and the numerical calculation 
comes from treating the energy injection rate as constant over 
9j and neglecting the contribution for 7 < 3. The above rough 
estimate is useful, which simply consists of the two factors, 
the ratio of the solid angle of the polar region to 4ir and the 
ratio of the energy in the effective injection range f\ < 7 < ri 
to the total range 1 < 7 < 7m- 

We confirm that the difference of the torus thickness tan 6q 
does not substantially change the value of the geometric effi- 
ciency. This is attributed to the fact that the energy injection is 
effective for large radii, for which the structure of the emission 
region is not relevant. The differences of the other parameter 
values, p,jM,ri, and tan(?i, affect the efficiency according to 
the simple formula Eq. ( fT2l . 

In fac t, 7m is typically much larger than 10 2 (see Sec- 
tion 12. Il l, which may significantly reduce the efficiency. To 
check this, we plot the factor fP -y~ p+l d~f / J^ M j~ p+l d~/ in 
Eq. ( flZb for f\ = 3, r-i = 10, and various 7m in Figure [5] It 
is found that this factor for 7m > 10 4 still can have a high 
value ~ 0.25 when the power-law index is p ~ 2.3-3. For 
7m = 10 6 and p = 2.5, the result of the numerical integration 
is 

-^- = 0.059 f n f a U, (13) 
Mc 2 

while the approximate estimate gives us Lj/(Mc 2 f n f a f t h) ~ 
0.04. Note that this corresponds to the maximum level of the 
efficiency in this system, which is obtained for the parameters 
p = 2.3-3,tanf?i = 1, and r2 = 10. However, we consider that 
such parameter values are within realistic ranges. 

The numerical calculation gives us the mass injection effi- 
ciency for 7m = 10 6 and p = 2.5 



M 



■■ 0.019 f„f a f th . 



(14) 
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This indicates that the injected protons will be isotropized and 



FIG. 6. — Results of the calculations of the energy injection efficiency eg = 
Lj/Mc 2 (filled circles) and of the mass injection efficiency em = Mj /M (open 
circles) as a function of Mg = M/(10 8 Mq) with the parameters fnfafth = 
3 X lfT 2 , tan O = 2, p = 2.5, -y M = 10 6 , r 2 = 10R S , and tan^ = 1. For M 8 > 

0. 3, eg oc My P+2 (represented by solid line) and em oc M g p+i (dashed line), 
consistent with the approximate formula Eq. H2\ . while for Mg < 0.1, ee 
and em both oc M 2 (solid and dashed lines), consistent with the approximate 
formula Eq. l !15t . The final Lorentz factor of the jet (without an additional 
energy input) is given by T = ee/em- 

have mean random Lorentz factor (7) = Lj/ (MjC 2 ) = 3 . 1 . If no 
additional energy is injected into the polar region, the fireball 
model predicts that the random kinetic energy is transferred to 
the bulk kinetic energy with terminal Lorentz factor r ~ 3 . 1 . 

4.2. Dependence on BH mass 

We also perform calculations for various values of BH 
mass. For the case of Mg = 10, the neutron emission region 
is set between p m = 10 and pm = 30. We find that many of 
the properties for this case are just a scale-up version of those 
for Ma = 1 discussed above. The energy is injected mainly 
for 7 ~ 30- 100 in the case of r2 = 100. The energy injection 
efficiency is approximately estimated by Eq. (I121 l, For p > 2, 
the energy in the range 7 ~ 30- 100 is smaller than that in 
7 ~ 3- 10, so that we have smaller energy injection in this 
case. Since the mean Lorentz factor of the injected protons 
{7) is larger, the mass injection efficiency is even smaller. 

We plot the calculated Lj/(Mc 2 ) and Mj/M for M 8 = 3 
and 10 with the same parameters as for Eqs. (TTjt and ([T4l 
and /,,/n/th = 3 x 10~ 2 in Figure [6] The approximate formula 
Eq. ( fT2l implies that Lj/Mc 2 oc r^ ,+2 oc M^ p+2 for given 9\ 
and 7m, which agrees with the numerical results. The mean 
Lorentz factor should linearly depend on the length scale, i.e., 

(7} oc M%, leading to Mj/M = Lj/((-f)Mc 2 ) oc Mj 1 ^ 1 , which 
also agrees with the numerical results. These scalings are ap- 
plicable for cases of r\ > 1, i.e., Mg > 0.3. 

For the case of Mg = 0.1, the system size is smaller than the 
case of Mg = 1, i.e., p m = 0.1, pm = 0.3, and ri = \ . In this 
case most of the energy is injected into a region at r, > r^, 
which is assumed to be significantly collimated, so that the 
injection efficiency should be smaller than the case of Mg = 

1 . Let us assume that the dilute polar region is cylindrical at 
rj > r2 for simplicity, and then integrate and ii[ p) /(jc 2 ) for 
the energy and mass injection, respectively, over 1 < 7 < 7m 
and the volume of the polar region at r, > 1 to deduce the 
efficiencies. The results (and those for the cases of Mg = 0.03 
and 0.2) are plotted in Figure [6] 

We can also derive an approximate formula of the energy 
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injection efficiency for Mg < 0.1. In the cylindrical region, 
we have the opening angle as a function of rj as 9 C ~ 9\r<i/rj. 
The geometric efficiency may be estimated as 

L, 1 



fnfaUMc 2 2Vp M r p+l d 1 



5 f°° ( 9 2 \ f 
d 1 J drjrj ( -±\ j dQjr^Sirj-^h 



9 1 r 2 J ! 



/r 7-^7 



4 f 7M 7 -p+irf 7 



:0.03r 2 2 , (r 2 < 1) 



(15) 



ie the final results are obtained for p = 2.5,7m = 10 6 , and 
?i = 1. For p ~ 2.3-3 and 7m > 10 , the efficiency ranges 



where 1 

tan 9\ = 1. For ~ 2.3-3 and 7m > 10 , the efficiency range 
between 0.02^ and 0.05^. This equation clearly indicates 
that the energy injection is smaller for smaller M%(= 0.1 r 2 ) 
in the case of r 2 < 1 . The mean Lorentz factor may be esti- 
mated as (7) = / 1 °°7~''~ 1 ^7/ 7 _/, ~ 2 <i7 ~ (p+l)/p. Those 
approximate formulae agree with the numerical results of 
Lj/Mc 2 and Mj/M. 

5. CONCLUSION AND DISCUSSION 

The luminosity of the relativistic neutrons from geo- 
metrically thick AGN accretion flow is estimated as L„ = 
fnfafthMc 2 , and f„f a f t h can be as high as - 3 x 10" 2 , where 
the thermalization fraction of the accretion power / t h ;$ 0.3, 
the energy fraction of the accelerated protons in the thermal 
energy f a < 0.3, and the neutronization (and escape) fraction 
fn ;$ 0.3. Those neutrons escape isotropically, a fraction of 
which decay into protons at the dilute polar region and inject 
the energy and mass into the relativistic jet. We have calcu- 
lated this geometric fraction by setting a simple system con- 
sisting of the central BH, the accretion flow including the neu- 
tron emission region (at R s < p < 3R S ), and the polar region. 
The polar region is assumed to be a cone with an opening 
angle of 9\ ~ 0(1) at r\ = R s < rj < r 2 = 10/? s and changed 
to be a cylinder, at rj > r 2 as a rough approximation of the 
significant collimation. 

The results are plotted in Figure |6] for which the param- 
eters are chosen as tan#o = 2, r 2 = 10R S , tan#i = 1, p = 2.5, 
7m = 10 6 , and f„f a f t h = 3 x 10" 2 , and they are well approxi- 
mated by Eqs. (pi forMg > 0.3 and (B) forMg < 0.1. The 
efficiencies Lj/Mc 2 and Mj/M do not significantly depend on 
the torus thickness characterized by tan U , because the energy 
injection is effective for large radii, for which the structure of 
the emission region is not relevant. The efficiencies are larger 
for the larger polar region, i.e., larger r 2 and/or 9\, according 
to Eqs. ([T2l and (Q3). They are also weak functions of p and 
7m for 2.3 < p < 3 and 7m > 10 4 , which are considered to be 
realistic ranges. 

For M ~ 10 8 M Q , this process can produce a relativistic 
jet with L; ~ 2 x 10" 3 Mc 2 and Mj~6x 10" 4 M, leading 



to the final Lorentz factor r ~ 3 (if no other types of ener- 
gies are injected). Most of the relativistic neutron luminosity 
L n ~ 3 x 10~ 2 Mc 2 is injected outside the polar region, which 
may contribute to the dense non-relativistic disk wind energy. 
The existence of su ch energetic disk winds is n ot incompatible 
with observations (Tombesi et al. 201 1, 2012). The main part 
of the accretion power Mc 2 is carried by protons and magnetic 
fields, which may also contribute to the relativistic jet and the 
non-relativistic disk wind, or just be swallowed by the central 
BH. 

The observations suggest that the luminosities of AGN jets 
are broadly distributed over L „■ tnt < Mc 2 (e.gJFernandes et al.1 



120111 iPunsly & Zhangll20Tll iGhisellird et alJl2010)) . Ourre 
suits imply that jets with Lj - Xot < 2 x 10~ 3 Mc 2 may be pro- 
duced by the neutron decay process by itself. In this model, 
larger BH masses are associated with larger final Lorentz fac- 
tors of the jets forM > 1O 8 M . 

For AGN jets with higher luminosities, the neutron decays 
still can be a dominant mass loading process, whereas they 
can have additional energy inputs, such as Poynting flux. The 
mass loading My ~6x 10" 4 M leads to e.g., T ~ 50 for L iitot ~ 

3 x 10" 2 Mc 2 . If the additional energy input scales as L, / tot 0^ 
M, jets with larger T are associated with smaller M for M < 
10 8 Mq and with larger M for M > 10 8 M Q . 

We have simplified the configuration of the BH accretion 
system, approximated the amount of the created protons as a 
delta function of the traversing distance of neutrons (Eq. |9), 
and assumed the Euclidean space, to obtain orders of magni- 
tudes of the geometric efficiencies. This study has provided 
the useful approximate formulae Eqs. (fL2t and (1151 . More 
sophisticated formulations and calculations, taking into ac- 
count the collimating shape of the dilute polar region, radial 
dependence of the neutron emissivity, detailed process of the 
neutron decay, and geometry around the Kerr BH, are worth 
investigating in separate papers. 

While we have focused on the geometrically thick disks 
in this paper, the proton acceleration and the relativis- 
tic neutron production may be effective also in hot coro- 
nae above the geometrically thi n disks (e.g.. iDruryl 120121 
Ide Gouveia Dal Pino et all 120101: I Vievro & Romero! BO 1 21) . 
Similar scalings of the energy and mass injections to the out- 
flows are expected. In any cases, it is important to discuss the 
mechanism for regulating r ~ 10- 100 of AGN jets and the 
role of the neutron component. 
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